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A quantum system weakly interacting with a fast environment usually undergoes a relaxation 
with complex frequencies whose imaginary parts are damping rates quadratic in the coupling to the 
environment, in accord with Fermi's "Golden Rule". We show for various models (spin damped 
by harmonic-oscillator or random-matrix baths, quantum diffusion, quantum Brownian motion) 
that upon increasing the coupling up to a critical value still small enough to allow for weak-coupling 
Markovian master equations, a new relaxation regime can occur. In that regime, complex frequencies 
lose their real parts such that the process becomes overdamped. Our results call into question the 
standard belief that overdamping is exclusively a strong coupling feature. 



PACS numbers: 05.30.-d; 03.65.Yz; 76.20.+q. 

I. INTRODUCTION 

The dynamics of an isolated and finite quantum sys- 
tem consists of a reversible superposition of oscillations 
with (real) Bohr frequencies u>s. In order to understand 
the irreversible processes occurring in finite quantum sys- 
tems, such as relaxation to equilibrium or decoherence, 
one needs to take into account the interaction between 
the system and its environment. The weak-interaction 
limit together with the Markovian approximation already 
allow a good understanding of such irreversible processes 
and has some universal features. The generator of the 
evolution of the (reduced) density matrix of the sys- 
tem obtained by second-order perturbation theory (often 
called the Redfieldian) is not an anti-Hermitian genera- 
tor any more. Its eigenvalues r + ifi acquire a real part T 
describing irreversible decay to equilibrium. The imagi- 
nary parts of the eigenvalues are shifted Bohr frequencies 
£1 = ujs — Su>. The two shifts T and 8u>, normally increase 
(quadratically) as the strength of the coupling grows. 

We here propose to show that Markovian perturba- 
tive master equations such as the Redfield equation 
S S S B B allow for more than just describing 
the well known normal damping just mentioned. When 
the coupling strength is increased, it can happen at a 
critical value that a shifted frequency f2 vanishes and 
for yet stronger coupling goes imaginary. The pertinent 
eigenvalues r — are real and, interestingly, decrease 
with growing coupling. The resulting relaxation is non- 
oscillatory, i.e. overdamped. The principle purpose of 
this paper is to show that contrary to common belief the 
transition to overdamping is still compatible with pertur- 
bative treatment. In brief, overdamping can be a weak- 
coupling effect. 

All models to be studied here have Hamiltonians like 

H = Hs + Hb + SB, (1) 



where Hs and Hb respectively generate the free motion 
of the system and the environment (bath) while the in- 
teraction involves respective coupling agents S and B. 

It may be well to emphasize that the so-called rotating- 
wave approximation 0, , extremely useful as it may be 
for very weak damping, in particular in quantum optics, 
is definitely not allowable for strong damping and over- 
damping. Indeed, the rotating-wave approximation is 
based on the assumption that the Bohr frequencies of the 
system are very large compared to the system damping 
rate such that all "anti-resonant" terms can be time aver- 
aged out when writing the master equation in the interac- 
tion picture. But overdamping occurs precisely when the 
Bohr frequencies of the system become of the order of or 
smaller then the system damping rate. In a recent study 
of low-quality resonators , the rotating-wave approxi- 
mation was shown to be still affordable for overlapping 
resonances. But the Hamiltonians to be employed in the 
present paper must retain the "anti-resonant" terms that 
the rotating-wave approximation would suppress. 

A word on physical contexts where overdamping shows 
up is in order. One such is diffusion, a topic to be dealt 
with below (Section 1111(1 . Another one is temporal fluc- 
tuations in critical phenomena, described by time depen- 
dent Ginzburg-Landau equations without inertial terms 
[13; Ref. El describes a derivation of such a Ginzburg- 
Landau equation from an underlying unitary evolution 
of a "larger" system. 

The plan of the paper is as follows: In section [HI we 
solve the Redfield master equation for a two-level sys- 
tem interacting with a general environment. When the 
environment is made of harmonic oscillators (spin-boson 
model) , we show in subsection III Bl that the transition 
from normal damping to overdamping occurs at a crit- 
ical value of the coupling which can be made arbitrar- 
ily small and therefore accessible to perturbation theory. 
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For environment operators Hb and B modeled by ran- 
dom matrices from the so-called Gaussian orthogonal en- 
semble (spin-GORM model), we show in subsection III CI 
that weak-coupling overdamping is compatible with the 
exact dynamics computed numerically. In section lTTTl wc 
show that the transition from a non-diffusive to a dif- 
fusive regime, recently identified for a particle traveling 
in a spatially extended system while interacting with an 
environment, corresponds in fact to a transition from nor- 
mal damping to overdamping; that transition will turn 
out amenable to perturbative analysis. Finally, in sec- 
tion [W] we study the transition from normal damping to 
overdamping for a central harmonic oscillator interacting 
with a large collection of harmonic oscillators (quantum 
Brownian motion). We show that overdamping again al- 
lows for perturbative treatment, by comparison with the 
exact results known for this model. Conclusion are drawn 
in section Ivl 



II. DAMPED SPIN 

A. Hamiltonian and Markovian master equation 

Any two-level system has the Pauli matrices <j x ,a y , a z 
(together with unity) as a complete set of observables. 
If such a "spin" interacts with a general environment we 
may choose the Hamiltonian as 



H = ^a z + H B + a x B . 



(2) 



Inasmuch as the interaction a x B does not commute with 
the Hamiltonians for the uncoupled spin and bath, it 
allows for transitions between the unperturbed energy 
levels. Denoting the means of the spin observables by 

x(t) = Trp(t)a x , y(t) = Trp(t)a y , z(t) = Trp(t)a y (3) 

we write the Redfield equation as 0, 

z(t) = 2T(z{oo) - z{t)) (4) 
x(t) = -u y(t) 

y(t) = { ) x(t)-2Ty(t), 

with the time dependent damping rate T(t) and fre- 
quency Q(t) and the stationary inversion z(oo) 



a(t) = (B(t)B(0)) of the bath coupling agent B (For 
definition and properties see appendix ^J. 

The Markovian approximation consists in taking the 
upper bounds of the time integrals in © to infinity, such 
that the damping constant and frequency become time 
independent, r(oo) = T, Sl(oo) = SI . That approxima- 
tion is legitimate when the spin dynamics characterized 
by the rates wo,^,r is much slower than the decay of 
the bath correlation function a(t) and requires that we 
restrict the further discussion to times much larger than 
the bath correlation time. We may then rewrite as 



T = ^(a(u! ) + a(-u )) 



(6) 



si 2 + r 2 = ujI + — u 2 l dujv^r- „ 



6>(lu) 



z(oo) 



a(-u ) - a(uj ) 
a(-uj ) + a(uj ) 



with a(ui) the Fourier transform of a(t). The solutions 
of equations (@J in the Markovian limit read 



-2Tt 



z(t) = z(oo) + (z(0) - z(oo))e 

= x(0)r- g (0H sin(m)e -r t 



(7) 



si 

+x(0) cos (Sit) e 



-rt 



M x(o)((si 2 + T 2 )/u )~ y (o)T . _ rt 

y(t) = sm (Sit) e 



si 

-y(0) cos (Sit) e" 



-Ft 



The reduced density matrix p = | + xb x + ya y + za z can 
thus be written as a superposition of four modes, 



p(t) = Y, c ^)p^ Ht - 



(8) 



For normal damping, si = 0, S2 = — 2r, S3 = —T + 1SI 
and S4 = — r — 1SI. Overdamping occurs when 



SI 2 < 0, 



(9) 



T(t) 



and then the rates of JSJ) are given by Si = 0, s 2 = — 2T, 

s 3 = -r + ini and s 4 = -r - ini. 



B. The spin-boson model 

Taking the bath as a collection of harmonic oscillators 
[Hlll2j| we have for its free Hamiltonian and coupling agent 



f 1 , , , „, % ,_, o. iz wc iiavc iui us lieu nanmiumiui aiiu tun miiiy atan 

-2 / dr cos(w r) C(t) (5) 

n Jo N N 

2^ 4 f\-t\r>t\ H B = -Y(Pn+^lQl), B = Ye n Q n . (10) 

h Jq n—l n—1 



n(t) 2 + r(t) 2 = lu, 



r(t)z(oo) = J dr sm(uj T) D(t) . 

Properties of the bath are represented by the func- 
tions C(t) and D(t), respectively the real and imagi- 
nary parts of the equilibrium autocorrelation function 



We assume a quasi-continuum of bath frequencies u> n , 
employ a spectral function ytu)) = J2 n e n^K ~~ w )> anc ^ 
adopt Ullersma's choice [see [T^J and Appendix Ia]. 

2 na 2 w 2 



7(w) 



(11) 
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where a is the decay rate of the autocorrelator of the 
bath coupling agent and k an overall coupling strength. 
Thus equipped we can evaluate the rates in {£J. In the 
limits of high temperature, i.e. Phwo = hico/hsT <C 1, 
we get 



This is in contrast to the normal-damping case, accessible 
from the above by replacing \J — 1 — > +i, where the two 
slowest rates |Re[s2]| and |Re[sa]| increase as the coupling 
becomes stronger. 



r ^° 2- 



na 



PH 2 (a 2 +uj 2 ) 



KCXLl) 2 



(3h 2 



r z(oo) 



fi(cr 



/3/l 2 (a 2 +c^) 



Ka 2 wo 



wg (13) 



(14) 



here the limit u>o/a — > has been taken to remain consis- 
tent with the Markovian approximation; note that in the 
present section k has the dimension of an action, such 
that r is a rate. 

The critical value of k at which overdamping occurs 
is now found with the help of Eq. 10 by subtracting 
Eq.JTSJ to the power two to Eq. (JT3J). We find 



13^0 h 2 f3cu 



(15) 



If k > k c , and if k c is small enough to be treated by 
perturbation theory we have a selfconsistent theory of 
overdamping. Clearly, high temperatures are favorable 
for that theory to apply since the pertinent k c is sup- 
pressed by the factor PHujq <C 1. 

One might fear that our way of obtaining k c is not 
completely consistent if solely restricted to second-order 
perturbation theory because T 2 is of order k 2 while f2 2 + 
r 2 is of order k and does not include the k 2 corrections. 
That fear would be eased by the following argument. If 
we were to add 0(k 2 ) corrections to the right-hand sides 
of Eqs. |12ll3p . the results (|T3|l for the critical coupling 
would be generalized to series in powers of the lea ding 
terms displayed in JTBJ. (In fact, Jang et al. Ref.[l4| 
found SI 2 + T 2 — lo 2 (1 + An 2 /h 2 ); by recalculating k c , we 
again find Eq. lfT5)l if (3HlJq «1) 

Another look at the high-temperature rates reveals an 
interesting feature of overdamping. We have from JSJ 



si 

S2 

s 4 



= 



(16) 



An 

2k 2k 
2k 2k 

'Wp ~ Wp 



IK ZK / K c 

Wp + WV 1 ~ v~ 
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Most remarkably, the slowest relaxation rate of the spin, 
| Re [33] |, decreases when the coupling to the environment 
increases. For strong overdamping, k c /k -C 1, we even 
have 



S3 



4k 



(17) 



C. The spin-GORM model 



We retain the overall Hamiltonian J5J but modify 
the environment so as to let the free-bath Hamiltonian 
Hb and the coupling agent B be represented by ran- 
dom matrices from the Gaussian orthogonal ensemble 
(GOEV The resulting spin-GORM model was studied in 
Refs. [a, E|- We use the results of that work; in particu- 
lar, we adopt a unit of time that makes the Hamiltonian 
dimensionless and the bath correlation time of order h 
(see Ea. (fT§|) below). Specifically, we write 



H B = 



X 



B = 77- 



X 



(18) 



here X and X are random y x T GOE matrices with 
mean zero. Their non-diagonal (resp. diagonal) elements 
have standard deviation ctjvd = 1 (resp. <jd = V2)- The 
parameter 77 serves as a coupling strength. 

To study this model it is convenient to assume that the 
environment is initially in a microcanonical distribution 
with the (dimensionless) energy e. The autocorrelator of 
the bath coupling agent then reads 



a(e,t) N ^°° ^W^))^ 



(19) 



and has the Fourier transform 



N- 



rfh 
~2T 



(e + ho) 2 



(20) 



It may be well to note that we here meet Wigner's semi- 
circle law for the mean level density of the GOE. 

The general rates of the Markovian Redfield equation 
given in Eq. © can be evaluated and read 



r(e) 



2h 



(e~hu ) 2 + 



(21) 
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and 

n(e) 2 + T(e) 2 = ujl+iful 
;2 J( £ + huo) 2 -k 



(22) 



■ arctan 



2 J(e + hw ) 2 - I 
—ujq- 



arctan 



h 7T 



■ arctan 



-w 



■ arctan 



(e- 


I- hwo) + 


i 

2 


\A £ 


+ hu ) 2 


1 

4 




\- hwo) — 


_L 

■> 


\A e 
V 


+ huo) 2 


_ 1 

A 


(e- 


- hwo) + 


1 
2 


\A £ 


- f^o) 2 


1 

4 


(e- 


- Hujq) - 


1 
2 


\A e 


- ftw ) 2 


1 
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If f2 2 < 0, we have overdamping. To discuss that case, 
we momentarily set ft 2 = A — B where A = £1 2 + T 2 and 
is given by (|2*T)l and B = T 2 by (|2*27> . When 77 is large we 
could have overdamping because B > A, but then per- 
turbation theory may fail and our approach lose selfcon- 
sistcncy. However, since all terms in A (but none in B) 
carry explicit factors loq or ujq , and since the other quan- 
tities containing u>q (i.e. and arctan(.)) are bounded 
away from zero in the limit luq — > 0, it is always possible 
to choose ojq sufficiently small such that A < B for small 
r). This is illustrated in Fig. ^ 




FIG. 1: 77 = 0.2, e = and H = 1. This figure shows that the 
condition for overdamping can be satisfied at weak coupling 
if coo is sufficiently small. This is still true for any generic 
choice of e. 

The dependence of the smallest rates on the coupling 
strength is similar as in the spin-boson model. The rates 
|Re[s3]| and jRefsJ] [see Eqs. JjjJ], grow with the cou- 
pling constant r\ in the normal-damping regime, have a 
cusp at the transition, and then decay into the regime of 
overdamping, as illustrated in Fig. |2 

We have numerically solved the exact dynamics in or- 
der to verify that the perturbative equation predicts the 
correct dynamics for normal damping as well as for over- 
damping. The agreement is excellent as illustrated in 
Fig. |21 We can conclude that the spin-GORM model 
allows for overdamping at weak coupling. 





2 0.25 0.3 



FIG. 2: loq — 0.01, e = and h = 1. The upper figure shows 
that for a fixed and small value of wo their exist a critical 
and small value of the coupling r\ c ~ 0.14 above which over- 
damping occurs. The lower figure illustrates the qualitative 
change of the coupling dependence of the slowest relaxation 
rates when going from the normal damping regime to the 
overdamped regime. 



III. DIFFUSION MODEL 

We now consider a particle moving on one dimensional 
closed loop while interacting with an environment. The 
pertinent dynamics has been studied recently in Refs. 
H [H3- 0| by using the Redfield equation. A transition 
from nondiffusive to diffusive relaxation has been identi- 
fied. We shall here use the results of this study to show 
that the transition mentioned in fact is one from normal 
damping to overdamping. 

The Hamiltonian of the loop constituting the subsys- 
tem is represented by an N x N matrix 



H S = 



E 
-A 







\-A 



-A 
E 
-A 





-A 

E 






-A 



-A 





Eq 

-A 








-A 

Eq 

-A 



-A 







-A 
E Q / 



(23) 



NxN 



taken in the site basis \l), where I = 0, 1, . . . , N— 1 labels 
the N sites on the loop. The diagonal elements of Hs are 
the on-site energies of the particle while the offdiagonal 
elements generate hopping to neighboring sites. 
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FIG. 3: (Color online) Transition from normal damping 
to overdamping in the spin-GORM model. The full lines 
represent the exact dynamics of the three spin observable 
x(t),y(t),z(t) obtained numerically by diagonalizing the full 
Hamiltonian and the dashed lines represent the dynamics pre- 
dicted by the Redfield equation (second order perturbation 
theory). The two results give curves which are so close to 
each other that the dashed lines are almost invisible. The 
situation depicted here is the same as in Fig. |5|where r\ varies 
and uio = 0.01, e = and h = 1. As predicted by Redfield 
theory, the transition occurs at r\ c ~ 0.14. The initial con- 
dition is x(0) = V%/3,y(0) = 0,z(0) = 1/3. For the exact 
dynamics we have taken N = 3000 and a width of the initial 
energy shell 8e — 0.025. 



A weak interaction with an environment is described 
by the Redfield master equation. The correlation time of 
the environment is assumed much shorter than all char- 
acteristic time scales of the loop and therefore the corre- 
lation function of the environment can be modeled by 



%(t)=2QJ(t) 5 w 



(24) 



By using the Bloch theorem, the Redfield generator (con- 
taining A* 4 elements) can be simplified in N independent 
sectors (with N 2 elements) , corresponding each to a given 
value of the Bloch number q. For our finite loop period- 
icity yields q = u2tt/N, where n — 1, 2, . . . , N. By diag- 
onalizing a given sector we get N eigenvalues depending 
on q. The complete spectrum of the Redfield generator 
then consists of the A^ 2 eigenvalues obtained by varying 

As already mentioned, two relaxation regimes have 
been identified in this model. In the nondiffusivc regime 
all eigenvalues are complex with real parts of similar mag- 
nitude, proportional to the coupling constant Q, 



Re[s] 



(25) 



However, in a given sector (therefore at a given q) when 
the coupling term is increased beyond the value Q — 
27L4sin|, one of the A" eigenvalues separates from the 
other A^ — 1 ones. This eigenvalue is always real and is 
called the diffusive one. The diffusive branch is made of 
the diffusive eigenvalues of the different sectors. These 
eigenvalues have a smaller magnitude than the real parts 
of the nondiffusive eigenvalues. They therefore control 
the long time relaxation of the subsystem. The diffusive 
eigenvalues are given by 



2Q 2Q I (2HA . q 



(26) 



Lets define Q c = 2Msin-^ w 2-kKA/N. For Q < Q c 
no diffusive eigenvalues are present in the spectrum and 
the relaxation regime is nondiffusive [see Eq. Ij25(l ]. As 
soon as Q > Q c , at least two diffusive eigenvalues exist 
in the spectrum and the relaxation regime is called the 
diffusive regime. The two smallest diffusive eigenvalues 
controlling the long time scale relaxation are 



4tt 2 A 2 
' QN 2 



(27) 



Notice that the perturbative approach is consistent, be- 
cause Q c can be made as small as desired by choosing 
A/N small. 

It is already clear at this point that the nondiffu- 
sive (resp. diffusive) regime implies normal damping 
(resp. overdamping). Indeed, as for normal damping 
(resp. overdamping), the smallest relaxation rates in- 
crease (resp. decrease) with growing coupling in the non- 
diffusive (resp. diffusing) regime. Furthermore, as in the 
normal damping (resp. overdamping) regime, the small 
Redfield eigenvalues are complex (resp. real) in the non- 
diffusive (resp. diffusing) regime. We can make that asso- 
ciation even clearer if we assume the environment made of 
harmonic oscillators which we model by using Ullersma's 
spectral density [see appendix 10] . In this case, we find 
that at high temperature, the zero-frequency limit of the 
Fourier transform of the environment correlation func- 
tion is given by 



r ~i ^ v J H K 
lim a{u) = lim = — - 

w-*0 UJ TTp 



(28) 



Since our instantaneous-decay assumption l|24|) implies 
a(0) = - (29) 

we conclude 

Q = J (30) 
and thus find the diffusive-branch eigenvalues 



2« 2k / . .o. B xo 



(31) 



BA 2 2 
1 
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The similarity between these diffusive eigenvalues and 
the smallest eigenvalue of the spin-boson model in the 
overdamping regime [S3 in Eq. I|16(l ] is obvious, as is 
similarity between the real part of the small nondiffu- 
sive eigenvalues [ (12511 with 1|30[) ] and the real parts of the 
small eigenvalues of the spin-boson model in the normal 
damping regime [S3 and S4 in (|16Jl ]. 



B. Exact treatment 

The Hamiltonian (|32ll generates the Heisenberg equa- 
tions of motion 



P(t) 



N 2 



N 



+ E772 )Q(t)~J2 £ Mt) 



IV. QUANTUM BROWNIAN MOTION (QBM) 



A. Hamiltonian 

In this section we study overdamping in an exactly 
solvable model of Brownian motion. The model is made 
of a central harmonic oscillator interacting with an en- 
vironment which itself is a collection of harmonic oscil- 
lators (see, e.g. HEIEHI; 

these references will lead 
the reader to earlier work). The exact solution proves 
extremely valuable for our endeavor since it will be seen 
to yield, in the Markovian limit, precicely the same con- 
dition for overdamping as the perturbative treatment. 

We write the total QBM Hamiltonian as [2(| 



1 1 N 

H= 1 -(P 2 + c 2 Q 2 ) + i]T(P„ 2 + c 2 (Q„ - ^Q) 2 )(32) 

n=l n 



and thus have the system and bath parts 



1 N 2 



1 N 



n=l 



(33) 
(34) 



The coupling agents of system and bath read 



JY 



S — Qo, ; B — — e n Qn ■ 



(35) 



n=l 



The QBM Hamiltonian (221 is a sum of squares and 
thus manifestly positive. A not manifestly positive vari- 
ant of that Hamiltonian 0, 0] , discussed in Appendix 
IbI can be mapped onto the QBM Hamiltonian by a renor- 
malization of the bare frequency of the central oscillator. 
That observation allows us to use the exact results of 
Ref. f° r our present study of QBM. 



P n (t) = -wlQ n (t) - e n Q(t) 
$(t) = P(t) 

4i(t) = Pn(t). (36) 

The solution of H3fi|) can be written as 

N 

Q»(t) = ^(VW0v(0)+V(i)Pv(0)) (37) 

v=0 



Ut)- 



the indices \x and v step from to iV and Qq = Q, Pq = 
P. All A^ v {tys can be expressed in terms of the function 



N 



N 



-s-e^-i: 



(38) 



The zeros of g(z) yield the eigenfrequencies of Eqs. 

Assuming the bath frequencies to form a quasi- 
continuum we employ a spectral function 7(0;) — 
e n^( w n — u ) to replace the sum in l|3^|) by an integral, 



duj 7 „ . (39) 



We adopt an initial condition with statistical indepen- 
dence of central oscillator and bath, without restriction 
for the density operator p(0) of the central oscillator, 



Ptot(O) =p(0)- 



Z b 



(40) 



The time dependent density operator of the central os- 
cillator then obeys the exact master equation 



+y m {t)[Q,[p,p{t)] + ] 



(41) 



-jpd PP (t)[Q,[Q,p(t)]} 

+^d pq (t)[P, [Q,p(t)]}, 

with [•,•]+ the anticommutator. The drift and diffu- 
sion coefficients f pq (t), f pp (t),d P p(t),d pq (t) can be found 
in |l9|: they can all be expressed in terms of the quan- 
tity A(t) = Aoo(t). To get an explicit result for that 
amplitude we adopt Ullersma's spectral function, 



7(w) 



2 na uj 



2, ,2 



(42) 
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where a and k are the decay rate of the autocorrelator of 
the bath coupling agent and an overall coupling strength, 
both now of the dimension of a frequency. For that choice 
the amplitude in question takes the form 



Ait) = 



2T 



a 2 + n 2 + r 2 - 2Ar 
a 2 + n 2 - r 2 



e~ rt cos(Q*)) 



1 



,-rt 



A 2 



n 2 + r 2 - 2Ar n 



sin(nt) . (43) 



Here, the three rates (r, f2, A) control the exact dynam- 
ics; they are connected to the three model parameters 
(luq, k, a) by the characteristic equations 



A = a- 2I\ 



2 

LJ 



an = Q 

i2 , -n2 



(n 2 



T z + 2XT . 
n(X/a). 



(44) 



The coupling between central oscillator and bath is thus 
seen to shift the unperturbed frequency as ojo — > + iT 
and the unperturbed bath decay rate as a — > A. 

We should mention that the (diffusion) coefficients 
dp P (t) and d pq {t), in contrast to (the drift coefficients) 
f pq (t) and f pp (t), also depend on the temperature. 

As a final comment on the exact solution of the model 
we would like to add that, due to the initial condition 
(JlijJ, we have (P n ) — (Q n ) = and therefore get the 
mean displacement of the central oscillator from (|37ll as 



(Q(t)) = A(t)(Q(0)) + A(t)(P(0)) 



(45) 



Turning to the Markovian limit we assume that envi- 
ronment correlations decay fast relative to the time scales 
of the central oscillator. In technical terms, we require 



a,A> |r + ifi| 



(46) 



That Markovian limit does not imply weak coupling. The 
characteristic equations (|44|l now become 



A = a 



= n 2 + r 2 



2aT 



(47) 



lu 2 = (Q 2 +T 2 ) 
and entail the explicit results 

r = f, ^ 2 = ^-^, a = «. 

The master equation now reads, for times t 3> a -1 , 

m = -^\P 2 +Lo 2 Q 2 1 p{t)\ 



(48) 



(49) 



i 

h 

i 

1? 



r[Q,[P,p(t)}+] 

[Q, [Q,p(t)]] 

(iu 2 (Q 2 ) cq -(P 2 ) cq )[P, [Q,p(t)\] 



The exact expressions for the stationary second moments 
(Q 2 )eq and (P 2 ) cq are lengthly and can be found in [T^ |: 
they are completely characterized by the three rates (r, 
ft, A) and by the temperature. 

In the Markovian limit under study, the amplitude A(t) 
in 1)43(1 also simplifies to 



A{t) 



n 



e _it sin(fii), <>l/a. 



(50) 



We can now see that overdamping arises when Q becomes 
a pure imaginary number or equivalently when luq < T. 
The transition between normal damping and overdamp- 
ing occurs at £1 = 0, for the critical coupling 



• 



(51) 



That critical coupling will have to be compared with the 
one obtained perturbatively. 



C. Perturbative treatment 

In order to compare exact and perturbative results we 
now look at the Redfield master equation for the QBM 
Hamiltonian 18] 



kt) = --[P 2 + (n 2 p + T 2 p )Q 2 ,p(t)} 



2H 

- l -T p [Q,[P,p{t)]- 



(52) 



cq [ 



[Q,p(t)]] 



-l((fi 2 + r 2 )<Q 2 ) cq -(F 2 ' 



cq 



where 



dt 



w o 

(f7 2 + r 2 ) = Wo 2 + 



D(t), 



2 



[P, [Q,p(t)}} 

(53) 

dt cosujot D{t) , 



2rp(-f > )cq 



dt cos coot C(t) , 



(^ + r 2 )(Q 2 ) cq -(F 2 ) oq = / dt 

Jo 



sinujot 



C(t). 



The Markovian approximation consists in taking the up- 
per bounds of the time integrals of l|53l) to infinity and 
is justified when the free motion of the central oscillator 
(characterized by the frequency uiq) is much slower than 
the characteristic decay rate a of the correlation function 
of the environment (uiq /a — > 0). Again using Ullersma's 
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spectral function l|42[l we get the foregoing rates as 



Acknowledgments 



p 2(a 2 +uj 2 ) 2 
(Q 2 p + T 2 p ) = uj 2 + kol 



(54) 



or 

a 



a 2 + uj'I 



2 fkJo ^fihujQ fj^a 1 



(n eq = -fcoth , 



(n p + T 2 )(Q 2 ) cq -(p 2 ) 





2\ f 3 Zl 



(3{a 2 +uo 2 ) (3 



To be consistent with the Markovian assumption, all 
terms of order ujo/a or smaller should be disregarded. 
When using the lowest-order master equation H52[) we 
recover the mean displacement (Q{t)) of the rigorous 
treatment; in fact, we even get coinciding results for 
the non-perturbative and the perturbative rates in the 
Markovian limit u>o/a — > 0, i.e. r — T p = k/2 and 
fl = H p = luq — k 2 /4. In particular, therefore, the transi- 
tion to overdamping occurs at the same critical value of 
the coupling, given by Eq. (|51|l . We conclude that the 
overdamping regime in the QBM model in the Marko- 
vian limit can be described by second-order perturbation 
theory and therefore is a weak-coupling overdamping. It 
is worth mentioning that this result was anticipated by 
Cohen- Tannoudji in [2l|. 

We finally note that for strong overdamping the slowest 
decay rate of the QBM model reads 



(55) 



in obvious similarity to the corresponding limit for the 
other models studied above [see and JST}]. 



V. CONCLUSION 

For four different models, made of a system weakly 
interacting with its environment, we have studied the 
transition from normal damping to overdamping. Nor- 
mal damping has slowest relaxation rates that increase 
with growing coupling strength and is characterized by 
exponentially damped oscillations. In the overdamped 
regime the smallest relaxation rates decrease with grow- 
ing coupling and the dynamics displays non-oscillatory 
exponential decay. The critical value of the coupling at 
which the transition from normal damping to overdamp- 
ing occurs can often be made sufficiently small (by tuning 
model parameters) to be describable by weak-coupling 
master equations such as the Redfield equation. One way 
to make the critical coupling small is to decrease the bare 
frequencies of the system, but other parameters like the 
temperature of the system size can also enter the game. 
The compatibility of weak coupling and overdamping is 
counter to intuitive and widely spread expectations. 
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APPENDIX A: HARMONIC OSCILLATOR 
ENVIRONMENTS 

We briefly recall some properties of the equilibrium 
autocorrelator of the environment coupling agent B, 

a(t) = (B(t)B) — C(t) + \D(t) (Al) 
= Tr B ^e- l6Bt ' h Be lflBt ' h B , 



Pb 



(A2) 



The real and imaginary parts of a(t) obey C{t) = C(—t) 
and D(t) = —D{—t). Their Fourier tarnsforms (defined 
as a(w) = ^ dte iLOt a(t) = C{uj) +lD(w)) are related 
by the fluctuation-dissipation theorem 



where 



. . hu> /3huj 
^s(w) = — coth — 



(A3) 



(A4) 



is the thermal energy of an oscillation with frequency to. 
As a consequence, we can write our correlator as 

a(t) — J dujhJ{uj) ^coth^^ cos uit — i sin totj , (A5) 

thus introducing the spectral strength J(uj) of the envi- 
ronment often used in the literature, 



J(uS) = %b{u)), oj>0. 
h 



(A6) 



It is in fact customary to use that spectral strength only 
for positive frequencies; an extension to real frequencies 
could be to require J to be odd in u>. 
For an oscillator bath with 



1 N — , 

11 '•■ = 2 ^ {p " + ' ^ = E e « ( 

71 = 1 



N 



71=1 

the correlator becomes 



N 



c -fiH B 



71=1 

N 



Z B 



Qn(t)Q n {0) 



- — I coth — - — cos ui n t — l sm uj n t I (A7) 
2cu n V 2 / 

71=1 



^coth— — cos u;t — isinwt^ 



2u) 
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and has the Fourier transform 

aluj) = ni u coth- + 1). A8 

Au 2 

A comparison of the general form (| A5|) with the 
oscillator-bath form (|A7|) of the correlator a(t) shows 
that the two spectral strengths J(ui) and j(co) (which 
are both common currency) are related as 



JM = lM, w >o. 



Ullersma's choice 0] (also called Drude strength) 

2 na 2 u) 2 



w a + uj 



(A10) 



corresponds to an ohmic environment because at small 
frequencies J {to) ~ klo/t:. 

At high temperatures, the real part of the environment 
correlator is given by 



C(t) 



7(w) 



KO> 

duj - „ cosuit = — — e 
o Pu 2 P 



-a\t\ 



(All) 



The imaginary part of the environment correlation func- 
tion is independent of temperature and reads 



D(t) 



7(w) . 



2w 



sin tot 



hua 2 



-a\t\ 



sgn(f) . 



APPENDIX B: ULLERSMA'S HAMILTONIAN 

Ullersma 01 an d other authors work with a mod- 
ified Hamiltonian of the oscillator model, 

1 1 N N 

H = -{P 2 +tolQ 2 ) + -YyPl+ulQl)+QY,tnQn. 



The potential-energy part 



N N 
V(Q, {Qn}) = - (ui 2 Q 2 + ]T W 2 n Ql) +QY, tnQn 



n=l 



has a minimum of the potential created by the other har- 
monic oscillators on the central oscillator given by 



(A9) dV(Q ,{Q n }) 



dQ n 



lQ„=<5„(min) = w nOn(min) + e n Q = . 



The central oscillator thus "feels" the potential 



V(Q ,{Qn(min)}) = (^-^^Ql 

n=l 71 

Clearly, then, positivity is not manifest; rather, in order 
to have bound states, we have to impose the condition 



J 2 



N 9 

\ - ei 



^ ,,,2 



lo — Ka > . 



(Bl) 



Ullersma's Hamiltonian can be mapped onto the QBM 
Hamiltonian by renormalizing the frequency ojq as lu 2 — > 
i^>o + = w o + Ka - That mapping was extensively 

used above in transcribing the rigorous results of Ref. 0] 
to the dynamics generated by the QBR Hamiltonian. 

Needless to say, we could have based our study of 
the transition from normal damping to overdamping on 
Ullersma's model. Only one subtlety about that alterna- 
tive treatment is worth being mentioned here. To leading 
order in coo/a the critical value k c of the coupling turns 
out to coincide with the border K ma x = <jJq /a to positivity 
loss following from (|B1J| . 



71=1 



n=l 
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